//Evaluate magnitude response of continuous-time system 
hs=analpf(4,'cheb1',[.1 0],5)
fr=0:.1:30;
hf=freq(hs(2),hs(3),%i*fr);
hm=abs(hf);
plot(fr,hm);

hst = analogTrans(hs, 5, 'bp', [10,20])
hft = freq(hst(2), hst(3), %i*fr);
hmt = abs(hft);
plot(fr, hmt);


